Liquids in multi-orbital SU(A^) Magnets with Ultracold Alkaline Earth Atoms 
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In this work we study one family of hquid states of fc— orbital SU(A'^) spin systems, focusing 
on the case of A; = 2 which can be realized by ultracold alkaline earth atoms trapped in optical 
lattices, with N as large as 10. Five different algebraic liquid states with selectively coupled charge, 
spin and orbital quantum fluctuations are considered. The algebraic liquid states can be stabilized 
with large enough A'^, and the scaling dimension of physical order parameters is calculated using a 
systematic 1/A^ expansion. The phase transitions between these liquid states are also studied, and 
all the algebraic liquid states discussed in this work can be obtained from one "mother" state with 
SU(2) X U(l) gauge symmetry. 



I. INTRODUCTION 



Spin liquid state as an exotic quantum ground state 
of strongly correlated systems has been studied for 
decadesii^. Thanks to the active search for spin liquids 
in materials in the last few years3>''*, people are encour- 
aged to believe in the existence of spin liquids in nature. 
The stability of spin liquid usually relies on large num- 
ber of matter fields which suppress the continuous gauge 
field fluctuations. For instance, in the famous organic 
salts n — (ET)2Cu2(CN)3^i^, one of the proposed candi- 
date spin liquid involves a spinon fermi surface, where 
the finite density of states of matter field tends to sup- 
press the U(l) gauge field^>^. When the spinon fermi sea 
shrinks to a Dirac point, one needs to introduce large 
enough flavor number {Nf) of Dirac fermions to stabilize 
the spin liquid. However, large Nf is difficult to realize 
in SU(2) spin system, therefore one is motivated to look 
for systems with large spin symmetries. Tremendous the- 
oretical and numerical efforts were made on SU(Af) and 
Sp(7V) spin systems with large N^—. 

It was proposed that spin-3/2 cold atoms can realize 
Sp(4) symmetry without fine-tuningl^. Recently it has 
been discovered that an exact SU(iV) spin symmetry with 
N as large as 10 can be realized with alkaline earth cold 
atoms without fine-tuning any parameter—. Because the 
electrons carry zero total angular momentum, all the spin 
components belong to nuclear spins and hence the inter- 
action between atoms are totally independent of the spin 
components, i.e. the system has SU(A'') symmetry with 

= 25' -I- 1 for nuclear spin-S*. Therefore the alkaline 
earth cold atom plus optical lattice is a very promising 
system to realize the long-sought spin liquids. Besides 
the SU(A'') spins, there is another orbital degree of free- 
dom associated with the alkaline earth atoms, because 
both the ^5*0 and ^Pq orbital levels (denoted as g and e 
respectively) have SU(A^) spin symmetryiS. 

Most generally this system has symmetry SU(A/)s x 
U(l)c X U(l)o. U(l)c corresponds to the conservation of 
the total atom number i.e. the charge U(l) symmetry; 
U(l)o corresponds to the conservation of rte — Ug i.e. the 
orbital U(l) symmetry. We will tentatively assume the 
system has an extra orbital Z2 symmetry corresponding 



to switching e and g i.e. exp(i-|fT^), therefore we take 
the hopping amplitude of these two orbitals to be equal, 
also the two intraorbital Hubbard interactions are equal. 
Weak violation of this Z2 symmetry will be discussed 
in this paper, and we will show that it is irrelevant to 
the main physics discussed in this paper. Under these 
assumptions, after straightforward algebraic calculations 
the Hamiltonian in Eq. 2 of Ref^iS, can be rewritten as 



a 

= 1 ■ ■ ■ N , and a = e, g. 



H.c.+^U{f 



(1) 



Here 



E 



arn "'lam 

is the total number of the atoms on each site, T° = 
'^lam'^ap'^iP™. IS the pscudospin vector of orbital levels. 
U, J and Jz are simple linear combinations between U, 
Vex and V in Ref.— , which are from the s— wave scatter- 
ings between atoms. Since different orbital channels have 
different scattering lengths, J and Jz terms are allowed 
to exist because otherwise the system will have an un- 
physical SU(2A^) symmetry. Eq. [T]is the starting point 
of our study, and since all the fermionic alkaline atoms 
under study carry half integer nuclear spins, N will be 
taken to be even hereafter. 

In order to obtain more solid and quantitative results, 
we will keep both orbitals of the atoms at half-filling 



= and put this model on 



a honeycomb lattice with only nearest neighbor hop- 
ping. Therefore on top of the global symmetries dis- 
cussed before, there is another particle-hole symmetry 
with Cjam — >■ VjCjam' ^^'^ = 1 ^ud —1 with j belonging 
to sublattices A and B respectively. If t is the dominant 
energy scale of the Hamiltonian, the half-filled fermions 
on honeycomb lattice is a semimetal with two Dirac val- 
leys in the momentum space located at Q = (±4|^,0), 
and at low energy the band structure can be described 
by the Dirac Lagrangian 



L : 
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■^alf,d^llJa, (70, 71, 72) = (t^ 



(2) 



The 2x2 Dirac matrices r* are operating on the two sites 
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in each unit cell on the honeycomb lattice. The Dirac 
fermion has two Dirac points at the corners of the Bril- 
louin zone, therefore there are in total Nf = AN flavors 
of 2-coniponent Dirac fermions, with an enlarged 0(8A^) 
flavor symmetry at low energy, which will be manifest 
after we rewrite the Dirac fermions in terms of Majo- 
rana fermions. The short-range interactions between the 
Dirac fermions are irrelevant at the free Dirac fermion 
fixed point. 

In the following we will mostly be focusing on the 
Mott Insulator phase of Eq. [1] with U dominant. Mo- 
tivated by the spin liquid and weak Mott insulator 
K — (ET)2Cu2(CN)3^iii^, we want the system close to the 
Mott transition so that at short distance it still behaves 
like a semimetal, while at long distance the electron Ciam 
fractionalizes. With J = Jz = 0, the existence of a frac- 
tionalized phase close to the Mott transition on the hon- 
eycomb lattice was shown with a slave rotor calculation in 
Ref,^^, and the fractionalized spinon has the same mean- 
field band structure as the Dirac semimetal. In our sys- 
tem with nonzero J and Jz, various strongly correlated 
liquid states with coupled spin, charge and orbital fluc- 
tuations can be realized in different parameter regimes of 
Eq. [1] and all the liquid states can be obtained from the 
U(l) X SU(2) spin liquid that will be studied first. 

II. LIQUID STATES 
A. U(l) X SU(2) spin liquid, the mother state 

As the first example of liquid state, let us take both 
[/, J dominate t, while keeping = tentatively. 
In this case the symmetry of Eq. [T] is enhanced to 
SV{N), X U(l)c X SU(2)o. When U and J both dom- 
inate the kinetic energy, the system forbids charge fluc- 
tuations away from half- filling n = iV on each site, and 
also forbids orbital-triplet fluctuations, i.e. the low en- 
ergy subspace of the Hilbert space only contains orbital 
SU(2)o singlet. The Young tableau of the S\J{N) repre- 
sentation on each site has two columns with N/2 boxes 
each column, which is a large- A'^ generalization of SU(2) 
spin-1 (Fig. [T]2). The half- filling constraint on the low 
energy Hilbert space implies that one can do a local U(l) 
rotation on the fermions, which will be manifested by 
introducing a U(l) gauge field coupled to the charge 
degree of freedom as usual. The orbital-singlet constraint 
on each site implies that the local SU(2)o transformation 
will not change the physical state, and this local invari- 
ance can be described by a SU(2) gauge field coupled to 
the orbital indices of Ciam- 

More formally, one can introduce the bosonic U(l) 
slave rotor bi and SU(2) slave rotor, 2x2 matrix field 
hai3, as well as fermionic spinon fiam as foUowingi^: 

We will call b the chargeon and hais the triplon field. 
h is a group element of SU(2), with SU(2)lxSU(2)k 




FIG. 1: a, the Young tableau of the representation of N 
S\J{N) fermions on each site when orbital is constrained to 
be SU(2)o singlet, N has to be an even number, b, c and d, 
Feynman diagrams which contribute to the RG flow of the 
velocity anisotropy Eq. 1131 the solid square stands for the 
vertex cr^7fc9fe. 



transformation: h — )■ Aij^hAin. The SU(2)i symme- 
try is the physical SU(2)o symmetry of the orbitals, 
while the SU(2)i^ symmetry is a local SU(2) gauge sym- 
metry, which leaves the physical operator Ciam invari- 
ant with an accompanied SU(2) gauge transformation 
on fiam- f — > Mjj^f. The chargeon bi grants the 
spinon fiam a U(l) gauge symmetry as usual bi — )■ bie^^\ 
fiam — fiame~'^^\ and bi also carries the U(l)c charge 
i.e. bi will couple to the external electromagnetic field 
if the fermions Ciam were electrons. The properties of 
U(l) and SU(2) slave rotors were discussed in Reff^^ and 
Refi^° respectively, although the SU(2) slave rotors in 
Ref Jfl was engineered from a very different set-up. 

The U(l) and SU(2) gauge symmetry can be mani- 
fested by reformulating the hopping term of Eq. [1] using 
the decomposition of fermion operator Eq. [3l 

^ = E -tblbjfLhl,h,,pf,p + H.c. (4) 
<ij> 

And spinon fiam hops effectively in a band structure de- 
scribed by the following meanfield Hamiltonian 

^ = E ^m.j,ap)flfjp + H.c. 
<i,j> 

{U^j,af>) = {blb,hl^h,p^) (5) 

The value of {Uij^ais) should be solved self-consistently. 
If the self-consistent solution {Uij^ap) ~ ^ap, the U(l) 
and SU(2) symmetries are preserved by this meanfield 
solution. And the fluctuation on the meanfield solution 
is the gauge fields: Uij,ap ~ (C/jj,a/3)e"'°*^"^'=i ^'^"^ 
The dynamics of slave rotor b and hap are given by the 
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meanfield decompositions —tblbj{fj^hl^phjppfjp) and 

-th\aphjpfi {blbjfLfjn) respectively. 

In the Mott Insulator phase but close to the Mott 
transition, the spin model after second order t/U per- 
turbation will be very complicated. However, there is 
another self-consistent way of studying this system. Mo- 
tivated by the existence of the spinon fermi sea of weak 
Mott insulator n — (ET)2Cu2(CN)3, we assume here the 
weak Mott insulator phase is a phase in which the char- 
geon bi and triplon h are both gapped, and the fermionic 
spinon fiam hUs the same mean field band structure as 
the original fermions Ciam hi the semimetal phase with 
Nf = AN Savors of 2-component Dirac fermions at low 
energy (Eq. [5]), and then we can check the stability of 
this state. This spin liquid state corresponds to a mean 
field solution {b\bjh\^phjp(j) = USap, which preserves the 
U(l) and SU(2) gauge symmetry. After taking into ac- 
coimt of the U(l) and SU(2) gauge fluctuation, the low 
energy field theory of this spin liquid is described by the 
following 2-|-ld electro-weak theory like Lagrangian: 



2N 

0=1 



1=1 



(6) 



Here (70,71,72) = (t^, t^, — r^). ip is the low energy 

- 477- ■ 4tv 

mode of spinon /, and ipi = e^~^f, ip2 — e~^~^iTy f . 
Unlike Eq. [51 in Eq. [5] each Dirac fermion -0 is a four 
component fermion, because it contains both the Dirac 
indices and SU(2) gauge indices. 

The global symmetry of Eq. [5] is SU(2iV), which 
is a combined symmetry of SU(7V) spin symmetry and 
Dirac valley rotation, tp transforms nontrivially under 
translation, space refiection, rotation, time reversal, and 
particle-hole transformation as following: 



Tri : X ^ X + ^ ^ ^ ip, 

1 \/3 2^ 

Tr2 : x^x + -,y^y+—,il)'^e'-^>''tp, 



T : i -t, V ^ 7 VV', a 



^ — — A^ 



.,2..zj. „ /li , „ aI 



"y ■■ y^-y, V'^7 M^V', 02,^2 ^ -02, -A^, 



PH : Cjam Vjcl^^, ip -> j'^fi^tp^ap ^ -Up, 



A^ /l3 _^ _ 4I —43 



R-27r/3 : V' ^ e 



(7) 



/i" are three Pauli matrices that operate on the two 
Dirac valleys. Notice that the hexagons of the triangular 
lattice form a triangular lattice with three sublattices, 
and Pa.x is the refiection centered at sublattice a of the 



three sublattices. Vectors fi = (0,1), r2 = (^,— ^), 

ra = (— i). In the equation above, transforma- 
tions of gauge field components are not shown unless they 
transform nontrivially. R27r/3 is the hexagon centered ro- 
tation by 27r/3. Notice that time reversal transformation 
(T) always comes with a complex conjugate transforma- 
tion, and hence T only changes the sign of the SU(iV) 
as well as SU(2)o generators that are antisymmetric and 
purely imaginary, therefore the SU(A^) and SU(2)o Lie 
algebras are preserved. 

The gauge symmetry and global symmetry together 
guarantee that none of the apparently relevant perturba- 
tions like fermion bilinears exists in the Lagrangian Eq. 
El When N is large enough the Lagrangian in Eq. [HI is 
a conformal field theory (CFT). The ellipses in Eq. [51 
include all the gauge invariant four fermion interaction 
terms which break the SU(2A'^) global symmetry down 
to the symmetries of the microscopic Hamiltonian Eq. 
[H All these four fermion interactions are irrelevant for 
large enough N. This CFT fixed point is a pure spin liq- 
uid state because both the charge and orbital fluctuations 
are forbidden. The scaling dimension of gauge invariant 
physical order parameters at this CFT flxed point can be 
calculated using a systematic 1 /N expansion in a similar 
way as Kef.—<^^^, with the results: 

Ae.[V5^] = 2+^, A,4^r.t^]=2-^. (8) 

Here TJ^ is the generator of the SU(2A^) flavor symmetry. 
SU(2Ar) current operators ip^f^TJ^i^ gain no anomalous 
dimension from gauge fluctuations. The order parame- 
ters of many competing orders are classified as fermion 
bilinears of this spin liquid states. For instance, the three 
sublattice SU(Af) columnar valence bond solid order pa- 
rameter is tjj^'^'ip plus two other degenerate configurations 
after translation along x direction. The SU(A^) ferromag- 
netic and antiferromagnet order parameter are 'ipjoT'^ip 
and ^pfi^T'^ip, with a = 1- ■ ■ N"^ — 1. We can see that the 
VBS and the AF order parameters have the same scaling 
dimension, and it is smaller than the scaling dimension of 
FM order parameter based on the 1/iV expansion. When 
the four fermion interaction is strong enough there is a 
transition towards a phase characterized by one of the 
fermion bilinear order parameters. 

We took Jz = at the beginning of this section, but 
the algebraic spin liquid discussed here is stable against 
small Jz , because Jz will not introduce any gauge invari- 
ant relevant perturbation to the fleld theory Eq. [Bl For 
instance, all the fermion bilinears are ruled out by gauge 
symmetries and symmetries in Eq. [71 already. There- 
fore a small Jz only renormalizes four fermion terms to 
Lagrangian Eq. [51 More physically speaking, turning on 
small Jz will not allow any orbital triplon state in the low 
energy Hilbert space, therefore the U(l) x SU(2) gauge 
field formalism is still applicable. 

For the same reason as the previous paragraph, if 
we introduce a small perturbation that breaks the or- 
bital Z2 symmetry exchanging the two orbital levels. 
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no extra relevant gauge invariant perturbations on Eq. 
m can be induced. This is simply because that the 
spinon ijj does not carry any physical orbital charge, 
therefore a small Z2 symmetry breaking will not be 
reflected in the CFT. For instance, in the semimetal 
phase the Z2 symmetry breaking will lead to a veloc- 
ity anisotropy between two orbitals: ^tcla^cj. 
Written in terms of fractionalized quantities, this term 
reads J2<i.j> ^tbjbjfj^hl^^af^^h^^isf.f}, which breaks the 
global SU(2)o — SU(2)i symmetry but still preserves 
the SU(2) gauge symmetry. The linear order effect of 
the St term on the spinon band structure is propor- 
tional to {b\bjh\^^a^^^j^hii,p), and this expectation value 
is evaluated in the spin liquid state. Since the spin 
liquid state is invariant under the Z2 symmetry e*"2°' , 

i^i^j^lafj.'^fiiy^i'^l^) — 0, hence at the linear order the 
band structure of / is unchanged. In fact, the veloc- 
ity anisotropy X]<i j> ^^'^l^^'^j '^^^ lead to the following 
gauge invariant but "Lorentz symmetry" breaking cou- 
pling in addition to the field theory Eq. (6] 

3 3 ; 

m=l 1=1 

where fc is a; or y. Tr[ft,^cr^/icr™] is odd under the orbital 
Z2 transformation — > — cr^, or exp(i-|cr^). Therefore 
as long as the SU(2) slave rotor h remains gapped, this 
term only induces irrelevant term after integrating out 
the gapped h. However, as we will see in the next section, 
after the condensation of /i, an anisotropic velocity of the 
spin liquid will be induced, and we have to evaluate this 
anisotropy with RG equation. 

SU(2) gauge field has been introduced in SU(2) 
and more generally Sp(2A'^) spin systems with single 
orbita l^"'^^" — , but there the local SU(2) gauge symmetry 
is a transformation mixing particle and holes of spinous, 
and hence there is no extra U(l) gauge field as in Eq. El 
This particle-hole SU(2) gauge symmetry has no straight- 
forward generalization to larger nonabelian gauge sym- 
metries. In our case the SU(2) gauge field stems from 
the physical orbital degeneracy, and a straightforward 
generalization to SU(fc) gauge field with fc— orbitals can 
be made, as long as the Hamiltonian favors a total an- 
tisymmetric orbital state. In this case we can again de- 
compose Ciara aS Ciani = hhapfifini with h £ SU(fc). 

When k is large the SU(fc) gauge field tends to confine 
gauge charges, and controlled calculations are difficult. 
However, here SU(fc) gauge field fluctuation corresponds 
to the constraint of antisymmetric orbital state, which 
is analogous to large— S" of SU(2) spin system with anti- 
symmetric orbitals. Therefore the gauge confined phase 
could be a semiclassical spin ordered phase. 

The credibility of the U(l) x SU(fc) gauge field formal- 
ism can be tested in one dimension, where many results 
can be obtained exactly. For instance, one of the fixed 
points of fc— orbital SU(iV) spin chain is described by the 
Wess-Zumino (WZ) model of SU(A^) group at level fc2^. 
At the SU(iV)fe fixed point, the exact scaling dimension 



of the Neel order parameter is A = jy^jy+fc) ■ apply 
the U(l) X SU(fc) gauge field formalism to this spin chain, 
the first order 1 /N expansion gives the scaling dimension 
of Neel order A = 1 — which is consistent with the 
exact result. The equivalence between WZ theory and 
the constrained fermion was proved in Refj^. In one di- 
mensional spin chains, the WZ fixed point is usually not 
stable^'' with half-filhng, in the U(l) x SU(fc) gauge field 
formalism this instability is due to the relevant Umklapp 
four-fermion terms for arbitrarily large N . However, in 
2-t-ld all the four-fermion interactions are irrelevant with 
large enough N , therefore at the field theory level the spin 
liquid is more realistic in 2 -I- Id than 1-1- Id. 

Recently it was proposed that the most general ground 
state for SU(A^) Heisenberg magnet with fundamen- 
tal representation on each site is a gapped chiral spin 
liquids^. A chiral spin liquid can be obtained by sponta- 
neously developing time-reversal and refiection symmetry 
breaking fermion gap 'ijjip in the U(l)xSU(2) algebraic 
spin liquid, which will lead to the Chern-Simons topo- 
logical field theory for the gauge fields. 

The U(l) X SU(2) spin liquid state is very constrained, 
since both the half-filling constraint and SU(2)o singlet 
constraint are imposed on each site of the lattice. In 
the following we will study several other liquid states in 
the same system, which can be obtained from softening 
part of the constraints on the U(l) x SU(2) spin liquid 
state. Therefore the U(l) x SU(2) spin liquid state is the 
"mother" state of everything else in this paper. 

B. U(l) spin-orbital liquid 

Now let us take U large, while keeping J and Jz small. 
When U becomes dominant, the system forbids charge 
fiuctuations, but allows for coupled spin and orbital fluc- 
tuations. In this case we can just introduce chargeon 
h and spinon as Ciam = with a local U(l) 

gauge symmetry. Here the spinon is equivalent 

to J2 13 ^ia/s fii3m with fipm defined in the previous sec- 
tion. Therefore the new spinon does not carry any SU(2) 
gauge charge, but carries the physical SU(2)o charge. If 
the fermionic spinon fiam fills the same mean field band 
structure as the original fermions Ciam *-e. (blbj) is a 
constant, the low energy fleld theory of this spin-orbital 
liquid is described by the following 3D QED Lagrangian: 

4JV 

^qed = ^ V'oTmC'^M ~ 'i-a^)i>a H (10) 

a=l 

with global flavor symmetry SU(4A^). The existence of 
this phase has been shown with a U(l) slave rotor mean- 
field calculation in Ref J^. This type of Lagrangian has 
been studied quite extensively in the past, because sev- 
eral other spin liquid states also have the 3D QED as their 
low energy effective field theorj*iiiHi2i. It is well-known 
that when Nf = AN is larger than a critical number. 
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the 3D QED is a CFT^^. Since this CFT fixed point in- 
volves both spin and orbital degrees of freedom (but no 
charge fluctuation), we will call this CFT fixed point a 
U(l) spin-orbital liquid. 

In the well-known staggered flux state of SU(2) spin 
system, Nf = l^^'^^i^'* , while in the spin-orbital liquid 
states of alkaline earth atoms under study, Nf = AN can 
be as large as 40, therefore it is a much more promising 
system to realize this CFT. The first order 1/Nf expan- 
sion gives us the following results for the scaling dimen- 
sions: 

A,,S^] = 2 + ^H^, Ag.STgid^] = 2 - ^^(11) 

is the generator of the SU(A^/) flavor symmetry 
group. Again the S\J{Nf ) current "07/^7^^^?/' gains zero 
anomalous dimension. We can compare the U(l) gauge 
field formalism and expansion to the exact result of 
SU(2iV) chains in one dimension, and the 1/N expansion 
gives the exact result as WZ model at level k = 1. 

Now let us again introduce the Z2 symme- 
try breaking perturbation J2<i j> ^^'^I'^^'^j = 

J:<^,J>^MbJf^^^^>J'f^^^■ Unlike the U(l) x SU(2) 
spin liquid discussed in the previous section, since 
{b\bj) 7^ 0, now this Z2 symmetry breaking will in- 
troduce the following gauge invariant perturbation to 
the field theory Eq. [101 which cannot be absorbed by 
rescaling ip: 

5L = sTpa^^kidk - iak)^, (12) 

here k — x,y only includes the spatial coordinates. Phys- 
ically this term corresponds to the velocity difference be- 
tween the e and g orbitals, and it can be viewed as Eq. 
[9] after the condensation of SU(2) slave rotor h. We can 
evaluate the RG flow of this term at the order of 1/iV 
trough Feynman diagrams Fig. c and d, as was done 
in Ref.^^. And the result is that 

ds _ 16 
dlnl ~ IStt^TV*- ^ ' 

Therefore this perturbation is irrelevant under RG flow. 

Now if we gradually increase J in Eq. [1] finally the 
orbital triplons will be excluded from the low energy 
Hilbert space, and the U(l) x SU(2) spin liquid state 
discussed in the previous section becomes the candi- 
date ground state. The phase transition between the 
U(l)xSU(2) spin liquid and the U(l) spin-orbital liquid 
can be driven by condensing the triplon field hap, which 
can also be parametrized as /i = (j)oI + i(j)ia^ + i<j)2'y'^ + 
i(f>3(7^, is a real 0(4) vector, and u'' are Pauli matri- 
ces. Further we can define CP(1) field z — (zi,Z2)*, and 
zi = (j^o ~ i03j -32 = 02 — *0i- Now this phase transition 
can be described by the following Lagrangian: 

c = + \{d,^j2 + ^^\^\' + ■■■ (14) 

1=1 



with critical point r = 0. Lew is given by Eq. |6l After 
the condensation of z, all three SU(2) gauge field will 
be higgsed and gapped out, and the remnant gauge field 
is a^. Based on the definition of / in Eq. [31 the conden- 
sation of hai3 implies that /'^^ and / becomes equivalent 
after a global SU(2) rotation. This phase transition is 
beyond the Landau's theory, because neither side of the 
phase transition can be characterized by an order pa- 
rameter. For general SU(fc) gauge symmetry with fc > 2, 
condensation of matrix field hap always gaps out all com- 
ponents of nonabelian gauge fields. 

Since the fermion number in Cew is large, one can use a 
systematic 1 /N expansion to study the universality class 
of the transition Eq. 1141 and the large fermion flavor 
number will suppress the SU(2) gauge fluctuations. For 
instance, in the large-TV limit, we can view the SU(2) 
gauge field completely suppressed, then the transition 
described by Eq. [Til belongs to the 3d 0(4) universality 
class. Notice that other gauge invariant and symmetry 
allowed couplings between z and ip are at very high order, 
and hence are irrelevant at this transition. For instance, 
coupling |zp^/'7o'!/' violates the particle- hole symmetry. 
And the Z2 symmetry breaking term Eq. [5] is also irrel- 
evant at this transition due to its high scaling dimension 
from power-counting. 

The phase transition between the ordinary Dirac 
semimetal phase and the U(l) spin-orbital liquid phase 
can be driven by condensing the chargeon b in Eq. [3] 
described by Lagrangian 

C ^ Cqed + \{d^- ia^)b\^ + rfc|6|2 -|- . . . (15) 

Cqed is given by Eq. 1101 The condensation of chargeon b 
will higgs the U(l) gauge field a^, and release the charge 
fluctuation from the constrained Hilbert subspace. In the 
large- iV limit when the gauge field fluctuation is frozen 
by fermions, Eq. [15] is a 3d XY transition. The velocity 
anisotropy Eq. [12] is an irrelevant perturbation at this 
transition as well, because the fluctuation of b will not 
affect the RG equation Eq. [T3] at the order of 1/N. A 
similar metal and weak Mott insulator transition is stud- 
ied in Ref.'^""'^^, where the condensation of the chargeon 
rotor b kills the U(l) gauge field fluctuation, and drive 
the Mott insulator with spinon into an ordinary metal. 



C. SU(2) spin-charge liquid 

The next situation we will consider is to keep J large, 
and make and U small. In this case the system forbids 
triplon excitations, but charge excitations are allowed. 
We can introduce spinon f^"^^ as Ciam — ^iapfl'^lnJ and 
hajs is the same SU(2) slave rotor as introduced in Eq. 

('2) 

[3] while f^J^ is equivalent to bifiam- This spin-charge 
liquid state has meanfield solution {h\^phjpp) ~ 6a|i^ and 
at low energy can be described by Dirac fermions coupled 
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special for fc = 2, for QCD Lagrangian with SU(fc) gauge 
group with fc > 2, since it is impossible to form SU(fc) 
singlet cooper pair, the global symmetry is the apparent 
SU(2Ar) X U(l) symmetry. 

Again the SU(2) spin-charge liquid can be obtained 
from the U(l)xSU(2) spin hquid by "releasing" the 
charge degree of freedom, through condensing chargeon 
field bi in Eq. |31 The Lagrangian is similar to Eq. [121 



(18) 



After the condensation of 5, / and /'^-' are identical based 
on their definitions. The transition between the SU(2) 
spin-charge liquid and the ordinary semimetal phase can 
be described by condensing triplon Za from the SU(2) 
spin-charge liquid state, with Lagrangian similar to Eq. 



FIG. 2: Schematic phase diagrams with two tuning parame- 
ters. All the phase transitions denoted as A are described as 
a Higgs transition of U(l) slave rotor h, such as Eq. 1151 and 
Eq. 1181 phase transitions denoted as B are Higgs transition 
of SU(2) slave rotor hap, or CP(1) field z coupled with SU(2) 
gauge field, such as Eq. [T?] and Eq. 1191 phase transitions 
denoted as C are Higgs transition of SU(2) vector x coupled 
with SU(2) gauge field, for example Eq. 1231 The phase transi- 
tions D axe Higgs transition of spinon z\ coupled with gauge 
field j4p. Notice that in (o), Jz is weak, while in (c) Jz is 
strong and the constraint = is always imposed on each 
site. The multi-critical points in these phase diagrams are 
discussed in section HI. 



with only SU(2) gauge field with a QCD like Lagrangian 



2N 



3 

Cqcd = ^alt^idf, - Y] iAl^)i;a + 



1 = 1 



a 



(16) 



Since this state involves both spin and charge excitations, 
we will call it a spin-charge liquid. At first glance, the 
global symmetry in Eq. [16] is SU(27V) x U(l), but the 
true symmetry is actually Sp(4A^) D SU(2Af) x U(l), and 
this Sp(4A^) group is a subgroup of the 0(8A^) symmetry 
group of the Dirac fermions in the semimetal phase with- 
out coupling to any gauge field. The enlarged Sp(4A^) 
symmetry was discussed in Ref.^^ in the tt— flux state of 
Sp(2A^) magnets with the same fleld theory as Eq. [TOl 

The Sp(4A^) symmetry not only contains the explicit 
SU(2A^) flavor symmetry of Eq. [TBI but also involves the 
pairing channel of because now there is no U(l) gauge 
field, and the gauge singlet cooper pairs of V' are physical 
operators. When fc = 2, The physical order parameters 
have scaling dimensions 



l^qcA^A =2- 



32 



A,,d[v;r^v^] = 2- 



16 



(17) 



e SU(2A^), and there are other fermion pairing bilin- 
ears with the same scaling dimension as ^T"^^ due to the 
enlarged Sp(4A'^) symmetry. The enlarged symmetry is 



C ^Cgcd+ lid, 



(19) 



1=1 



of the CP(1) field, physical 
are identical after a global 



After the condensation 

fermion Ciam and spinon , 
SU(2) rotation. In the large- TV hmit, Eq. [THand Eq. [TH 
describe a 3d XY and 3d 0(4) transition respectively. A 
similar orbital Z2 breaking term is present in the field 
theory Eq. [16] and Eq. [TO] but again this term only leads 
to irrelevant effects. 



D. U(l) X U(l) spin-orbital liquid 

If J is small compared with t, while both U and are 
much larger, then although the charge fluctuation will 
still be forbidden, the Hamiltonian gives a green light 
to one component of the orbital triplet state: the state 
(|e,g) + \g,e))/V2 with = 0. Therefore there are 
two U(l) constraints on the system: Ue + ng — N, n,, — 
Ug — 0, therefore we need to introduce spinon which is 
invariant under both U(l) charge rotation and orbital 
rotation generated with . Therefore in the proximity 
of the semimetal phase, the most natural liquid state with 
these constraints on the honeycomb lattice is described 
by the Lagrangian with two U(l) gauge fields 

2^ _ 

C-qed2 = Y ^"-^t" ^^t' ^ - ^)'0a H (20) 

a=l ^ 

with flavor symmetry SU(27V)+ x SU(2A^)_ x Z-^. The two 
SU(2iV)± groups are generated by T/ = t£{\ ± ct^)/2 
respectively, and the Zi symmetry exchanges ±. The 
scaling dimensions of gauge invariant operators to the 
first order of 1/iV are 



If we turn on the Zi symmetry breaking perturbation on 
the lattice, as in the third section, the following term will 
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be induced in the field theory: 

5L = s-4)a^"fk{dk - iak)ii - '^iA^'i/j-fki'- (22) 

The scaling dimension of this perturbation can be calcu- 
lated in the same way as the third section, and it is still 
irrelevant according to the RG equation at the order of 
In the rest of this paper this Z2 symmetry breaking 
will not be mentioned unless it is relevant. 

Since the U(l) xU(l) spin-orbital liquid only allows one 
of the orbital triplet states, We can obtain the U(l) x 
U(l) spin-orbital liquid by higgsing two components of 
the SU(2) gauge field in the U(l) x SU(2) spin hquid 
discussed before. We already showed that condensing a 
fundamental spinor of the SU(2) gauge group will gap 
out all three components of the gauge fields, but if we 
just condense an adjoint vector of SU(2) gauge group, 
only two of the three components of the gauge field will 
be gapped. Therefore starting with the "mother" state 
U(l) X SU(2) spin liquid, the U(l) x U(l) spin-orbital 
liquid can be obtained by condensing SU(2) vector x = 
z^a'^z instead of z itself. This transition can be described 
by the field theory 

3 ^ 3 
C = Ce^+Y,- i^^X^ - ^^Mj^Xkf + ■■■ (23) 

tijk is the total antisymmetric tensor, and also the ad- 
joint representation of SU(2) gauge group: t^^ = icaij- 
Without loss of generality, we take x condense along the 
direction (0, 0, 1), then Aj^ and are gapped out, while 
remains gapless, which is the same as the U(l) x U(l) 
spin-orbital liquid. Notice that x is not a vector of the 
physical SU(2)l symmetry. To see this explicitly, we can 
rewrite x as 

X = z^az r^TT[h^a''hd]. (24) 

h is the SU(2) rotor introduced in Eq. [3] It is explicit in 
this equation that x is only invariant under the U(l) sub- 
group generated by T^, which is the physical symmetry 
of the system with finite Jz ■ Similarly, if we condense the 
SU(2) vector xi ^ Tr[/i'''<T^/i(T] from the mother state, we 
would obtain a state with constraint = on each site. 

In the large- limit the SU(2) gauge field is again sup- 
pressed by the fermions, and the transition Eq. [231 is a 
3d 0(3) transition. A similar phase transition was dis- 
cussed in a different context^. This field theory was also 
used as a trial unified theory of electro- weak interaction 
in particle physics, and the gapless A"^ was identified as 
the photon'^'^. However, nature chooses a different theory. 
For larger fc, condensing adjoint vector of SU(fc) gauge 
group always leaves some components of the gauge field 
gapless. 

E. U(l) spin-charge-orbital liquid 

Finally, if we only keep Jz large, while keeping U and 
J both small, the only constraint on the system is TJ = 



on each site. Then the candidate liquid state in this case 
is described by the following lagrangian 

2N 3 

Lqed3 = J2 ^^It^i^t^ - Y)V'a + • • • (25) 

a=l 

This state has spin, charge and two orbital states fluctua- 
tion, therefore following our convention this state will be 
called U(l) spin-charge-orbital liquid. This state can be 
obtained from the U(l)xU(l) spin-orbital liquid state by 
condensing chargeon 6. Moreover, this U(l) spin-charge 
liquid state can also be obtained from condensing SU(2) 
gauge vector x in the SU(2) spin-charge hquid Eq. [T6l 
which gaps out both A^^ and A^. In the large- limit, 
these two transitions belong to the id XY and 3(i 0(3) 
universality class respectively. 

We can also drive a direct transition between the U(l) 
spin-charge-orbital liquid and the semimetal phase, as 
long as we can gap out the U(l) gauge field in Eq. [511 
In the previous paragraph we mentioned that the U(l) 
spin-charge-orbital liquid state can be obtained from con- 
densing vector X in the SU(2) spin-charge liquid Eq. 1161 
After the condensation of x^ the degeneracy between the 
two CP(1) fields zi and Z2 is lifted, due to the gauge 
invariance coupling x ' (^z. Therefore zi and Z2 can 
condense separately, but not together. If one of Za con- 
denses, it will Higgs the gauge field A^, and drive the 
system into the ordinary semimetal phase. This transi- 
tion can be described by the field theory 

C - Cq,d3 + \{^^,~ ^4)^lP + r.ilzip + (26) 

In the large- limit this transition again belongs to the 
3d XY universality class. 

III. PHASE DIAGRAM AND MULTI-CRITICAL 
POINTS 

The phase diagram with two tuning parameters J and 
U and weak constant Jz is depicted in Fig. [2}i, with 
four different liquid phases. And there is a multi-critical 
point with both masses of Zq and b vanish. At this multi- 
critical point, the field theory reads 

L^L,^ + \{d, - za^)6p + \{d, - ^z4-)z|2 + . .(.27) 

In the large- limit z and b behave like a 3d 0(4) and 3c? 
XY transition respectively. On top of this field theory, 
the symmetry allows the interaction between b and Zq., 
such as |zp|6p. It is well-known that at the id 0(4) 
and XY transitions, the scaling dimensions A[|zp] > 
A[|5p] > 3/2^4, therefore this term |2p|5p is an irrel- 
evant perturbation in the field theory Eq. [271 Notice 
that in principle the velocity of b and z are different, and 
the velocities will flow under RG equation with flnite A^. 

Fig. [5]: is the phase diagram with two tuning parame- 
ters J and U and strong constant Jj, where the orbital 
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constraint = is always imposed. Again there is a 
multi-critical point with both masses of and b vanish. 
The field theory at this multi-critical point is 

3 ^ 3 

1=1 ^ j,k=l 

+ + (28) 

and it is clear that the coupling is irrelevant in the 

large-TV limit due to the fact that a[|xP] > A[\b\^] > 3/2. 

The multi-critical point in Fig. ^ with tuning param- 
eters Jz and U can be simply described by field theory 
\{dfj, — ia^)b\'^ + \ {dfj, — A^)zi\'^ , and it is stable against in- 
teractions between b and zi. However, the multi-critical 
point in Fig. [2li is no longer a simple combination be- 
tween Eq. [Ml Eq. [19] and Eq. [26l because the symmetry 
of the system allows the coupling x • z^az. The fate of 
this relevant perturbation is unclear at this point. 

IV. SUMMARY AND EXTENSIONS 

In summary, we studied five examples of liquid states 
motivated by the orbital flavor and large spin symmetry 
of alkaline earth cold atoms. The schematic phase dia- 
grams are depicted in Fig. [5] Experimentally the spin 
correlation calculated in this paper can in principle be 
measured using the momentum density distribution and 
noise correlation between atom spins proposed in Ref.'^^, 
after releasing the atoms from the trap. The VBS or- 
der which breaks the lattice translation symmetry also 
has algebraic correlation in the liquid states. The two 
atoms within one valence bond have stronger AF corre- 
lation J ~ t'^/U compared with other atoms, and hence 
tend to move closer to each other from the minima of the 
wells. This super-lattice structure can also be measured 
by the density correlation between atoms, which can be 
detected by the noise correlation^^. 

It would also be interesting to test the results of this 
work by numerically simulating model Eq. [l]as in Reffi^, 



since the system is fixed at half-filling, the sign prob- 
lem of ordinary interacting fermion system is no longer a 
concern. Analytically it is useful to pursue a slave rotor 
meanfield calculation like Ref J^. The liquid states dis- 
cussed in our paper is expected to occur in a finite region 
close to the Mott transition, and by tuning U still larger 
there can be a transition from our states into a state 
with background nonzero gauge flux through plaquette, 
or dimerized valence bond solids^^. In our case the inter- 
play between the U(l) and SU(2) slave rotors make the 
meanfield calculation more complicated, and more mean- 
field variational parameters need to be introduced. We 
will study this calculation in new future. 

In the current work, the universality class of all the 
phase transitions between different liquid states was only 
discussed in the large- Af limit. Since the number of bo- 
son field at this transition is not large, the 1/A^ expansion 
at the transition is actually nontrivial. Let us take the 
simplest quantum critical theory Eq. [T5| as an example. 
If there is no gauge field a^, the transition is a 3d XY 
transition, whose critical exponents can be obtained by 
summing over the Feynman diagrams to all orders of an 
e = 4 — d expansion. The 1 /N correction from the gauge 
field propagator will enter the Feynman diagram at ev- 
ery order of e expansion, therefore it is a nontrivial task 
trying to sum over all the diagrams at 1/A^ order. How- 
ever, if we generalize the boson number to large Ni,, then 
a systematic expansion of both 1/N and l/Ni, can be 
straightforwardly carried out, as was studied in Ref.^^. 

Our formalism can be applied to other multi-orbital 
magnets, including transition metal oxides with orbital 
degeneracy. We will explore this possibility in future. 



Acknowledgments 

The author thanks M. Hermele and S. Sachdev for 
helpful discussions. This work is supported by the So- 
ciety of Fellows and Milton Funds of Harvard University, 
and NSF grant DMR-0757145. 



1 P. W. Anderson, Mat. Res. Bull 8, 153 (1973). 
^ R W. Anderson, Science 235, 1196 (1987). 
^ Y. Shimizu, K. Miyagawa, K. Kanoda, M. Maesato, and 
G. Saito, Phys. Rev. Lett 91, 107001 (2003). 

* Y. Kurosaki, Y. Shimizu, K. Miyagawa, K. Kanoda, and 
G. Saito, Phys. Rev. Lett 95, 177001 (2005). 

^ O. I. Motrunich, Phys. Rev. B 72, 045105 (2005). 
® S.-S. Lee, P. A. Lee, and T. Senthil, Phys. Rev. Lett 98, 
067006 (2007). 

N. Read and S. Sachdev, Phys. Rev. B 42, 4568 (1990). 

* S. Sachdev and N. Read, Int. J. Mod. Phys. B 5, 219 
(1991). 

® R. Flint, M. Dzero, and P. Coleman, Nature Physics 4, 643 
(2008). 

° R. Flint and P. Coleman, Phys. Rev. B 79, 014424 (2009). 



" M. Hermele, T. Senthil, and M.P.A. Fisher, Phys. Rev. B. 

72, 104404 (2005). 
^2 Y. Ran and X.-G. Wen, Cond-mat/0609620 (2006). 

C. Xu and S. Sachdev, Phys. Rev. Lett 100, 137201 (2008). 

C. Xu, Phys. Rev. B 78, 054432 (2008). 

F. F. Assaad, Phys. Rev. B 71, 075103 (2005). 
^® A. Paramekanti and J. B. Marston, J. Phys. Cond. Matt. 

19, 125215 (2007). 

C. Wu, J. P. Hu, and S. C. Zhang, Phys. Rev. Lett 91, 
186402 (2003). 

^« A. V. Gorshkov, M. Hermele, V. Gurarie, C. Xu, P. S. 
Julienne, J. Ye, P. Zoller, E. Demler, M. D. Lukin, and 
A. M. Rey, arXiv:0905.2610 (2009). 

S.-S. Lee and P. A. Lee, Phys. Rev. Lett 95, 036403 (2005). 
^° M. Hermele, Phys. Rev. B. 76, 035125 (2007). 



9 



21 


W. Rantner and X. G. Wen, Phys. Rev. B. 66, 144501 


30 


T. Senthil, Phys. Rev. B 78, 045109 (2008). 




(2002). 


31 


D. Podolsky, A. Paramekanti, Y. B. Kim, and T. Senthil, 


22 


I. Affleck, Z. Zou, T. Hsu, and P. W. Anderson, Phys. Rev. 




Phys. Rev. Lett. 102, 186401 (2009). 




B 38, 745 (1988). 


32 


D. A. Pesin and L. Balents, arXiv:0907.2962 (2009). 


23 


E. Dagotto, E. Pradkin, and A. Moreo, Phys. Rev. B 38, 


33 


H. Georgi and S. L. Glashow, Phys. Rev. Lett 28, 1494 


24 


R2926 (1988). 


34 


(1972). 


X. G. Wen, Phys. Rev. B 65, 165113 (2002). 


P. Calabrese, P. Parruccini, A. Pelissetto, and E. Vicari, 


25 


I. Affleck, Nucl. Phys. B 265, 409 (1986). 




Phys. Rev. B 70, 174439 (2004). 


26 


S. G. Nacuhch and H. J. Schnitzer, Nucl. Phys. B 347, 687 


35 


E. Altman, E. Demlcr, and M. D. Lukin, Phys. Rev. A 70, 


27 


(1990). 


36 


013603 (2004). 


I. Affleck, Nucl. Phys. B 305, 582 (1988). 


The author is grateful to Eugene Demler for pointing this 


28 


M. Hermele, V. Gurarie, and A. M. Rey, Phys. Rev. Lett. 




out. 




103, 135301 (2009). 


37 


R. Kaul and S. Sachdev, Phys. Rev. B 77, 155105 (2008). 


29 


T. Appelquist and D. Nash, Phys. Rev. Lett 64, 721 








(1990). 







